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Abstract 

Catalysis of nucleon decay in white dwarfs is used to constrain the abun- 
dance of magnetic monopoles arising from Grand Unified Theories. Re- 
cent discoveries of the dimmest white dwarf ever observed, WD 1136-286 
with L = lO~ 4 ' 94 Z/0, place limits on the monopole flux. An abundance 
of monopoles greater than the new bound would heat this star to a lu- 
minosity higher than what is observed. The new bound is (F/cm ~ 2 s _1 
sr _1 ) (ox>/10~ 28 cm 2 ) < 1.3 x lCr 20 (-u/lCr 3 c) 2 , where v is the monopole 
velocity. The limit is improved by including the monopoles captured by 
the main-sequence progenitor of the white dwarf: (F/cm ~ 2 s _1 sr _1 ) 
(o-i;/10- 28 cm 2 ) < 3.5(26) x 1(T 21 for 10 17 (10 16 ) GeV monopoles. 

We also note that the dependence on monopole mass of flux bounds 
due to catalysis in neutron stars with main sequence accretion has previ- 
ously been calculated incorrectly (previously the bound has been stated as 
-F(iTi;/10 _28 cm 2 ) < 10 -28 cm ~ 2 s _1 sr _1 ). We show that the correct bounds 
are somewhat weaker for monopole mass other than 10 17 GeV. 
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I. INTRODUCTION 



The question of whether or not magnetic monopoles exist has intrigued theorists and 
experimentalists for a long time ffl. In 1974, t'Hooft |2| and Polyakov H independently 
showed that magnetic monopoles always appear as stable topological entities in any Grand 
Unified theory (GUT) that breaks down to electromagnetism. Hence, if Grand Unified 
theories are shown to be correct, monopoles of mass in the range 10 15 - 10 19 GeV should 
exist. Rubakov M and Callan [0] calculated that these monopoles catalyze nucleon decay 
with a cross section characteristic of strong interactions, av ~ 10~ 28 cm 2 . 

The abundance of these monopoles is an open question. The Kibble mechanism predicts 
roughly one monopole per horizon volume at the time of the Grand Unified phase transi- 
tion. However, this estimate provides a severe overabundance of the number of monopoles: 
monopoles overdose the Universe by many orders of magnitude. Instead an inflationary 
epoch U may reduce their density in the Universe. Then the present abundance is difficult 
to estimate. A clue for experimentalists about what monopole flux to expect can be provided 
by astrophysics. 

The Parker bound on the flux of monopoles was obtained by requiring survival of fiG 
magnetic fields observed in our Galaxy and gave F < 10~ 16 cm~ 2 sr _1 sec^ 1 . Subsequent 
improvements on this work include consideration of the monopole velocities |5J due to accel- 
eration by the galactic magnetic field. Another improvement is the extended Parker bound, 
which required survival of a smaller seed magnetic field in the early period of the Galaxy 
g: F < 1.2 x l(T 16 ( cm" 2 s" 1 sr" 1 . 

Another class of methods for determination of the monopole flux is based on the hy- 
pothesis that GUT monopoles give rise to the catalysis of nucleon decay. The basic idea 
is that monopoles traveling through the Galaxy lose enough energy to be captured in an 
object (e.q. white dwarfs, neutron stars, etc.) where they subsequently catalyze nucleon 
decay. The energy produced by the nucleon decay heats up the object and results in a flux 
of photons from the surface of the object. One can then compare this predicted luminosity 
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with what is actually observed. One must ensure that the monopoles would not make the 
object brighter than what is seen. The coolest star (or other object) seen provides the tight- 
est limit on the monopole flux. If there were more monopoles than allowed by the bound, 
then the dimmest star observed could not exist. 

Several authors have carried out this kind of analysis in neutron stars IIJ , nearby pulsar 



and white dwarfs. The strongest bound was obtained from consideration of the catalysis 
process in PSR 1929+10, an old pulsar [|ll|]. From this pulsar, the bound on the product 
of monopole flux times cross section for catalysis is (F/cm _2 sr _1 sec _1 )(cru/lCr 28 cm 2 ) < 
7 x 10~ 22 . If one includes the monopoles captured by the main sequence progenitor of the 



white dwarf, this bound becomes even tighter |12| 



(F/cm -2 sr _1 sec _1 )(cru/10 -28 cm 2 ) < 10 -28 . The consideration of monopole dynamics 
inside superconducting neutron-star cores leads to a bound 5 x lO -24 ^ 2 cm' 2 sr _1 s _1 fL3"f , 
where Tio is the age (in 10 10 years) of the pulsar's present magnetic field. 

As neutron stars are the densest astrophysical objects observed, they give rise to the 
tightest catalysis bounds. However, there is a certain amount of uncertainty due to the fact 
that the interiors of neutron stars are not well understood. For example, neutron stars can 
have very large magnetic fields ~ 10 12 G of unknown topology, and the motion of magnetic 
monopoles inside the neutron star would undoubtedly be affected by these magnetic fields. 
In addition neutron star interiors may contain pion condensates, again with uncertain effects 
on the monopoles. Because of the uncertainties with neutron star interiors, we turn to the 
next densest astrophysical objects in the Universe, white dwarfs. These stellar remnants are 
far better understood. The flux limits obtained from consideration of the catalysis process in 



white dwarfs are therefore important. Previously Freese [|1J] considered monopole catalyzed 
nucleon decay in white dwarfs. By comparing with the lowest luminosity white dwarf that 
had been seen at that time, she obtained a limit 

(F/cm- 2 s" 1 sr" 1 )((Tt;/10" 28 cm 2 ) < 2 x 1(T 18 . (1) 

The present work is motivated by new observational data of cool white dwarfs |15 



In particular, Bergeron, Ruiz, and Leggett found a white dwarf 1136-286 (ESO 439-26) 
with luminosity 1O _4 ' 94 L ; this is the dimmest white dwarf observed to date. We use the 
measured luminosities of old white dwarfs to constrain the radiation due to monopole- 
catalyzed nucleon decay and thus to obtain an upper limit to the average flux of monopoles 
in the Galaxy. Since a white dwarf with luminosity 1O _4,94 L is observed today, we know 
that the monopole-induced contribution to the white dwarf luminosity cannot exceed this 
value. These new data improve the limit on the monopole abundance due to catalysis in 
white dwarfs Jill by roughly two orders of magnitude. Of course, as dimmer white dwarfs 
are found, the bound will continue to get more restrictive. 

A monopole flux saturating this bound would keep the white dwarfs at luminosities 
at least this great and would lead to the prediction that no cooler white dwarfs will be 
found. As we will discuss, if it were indeed true that monopoles are keeping dwarfs hot, one 
would expect a different dependence of white dwarf luminosity on mass than expected in 
the standard model. 

We shall explicitly indicate the dependence of our results on various parameters. We will 
parametrize the properties of the white dwarf in terms of typical values from observations: 
for the mass, M = M o e 0.6M & , for the radius R = Rg9 x 10 8 cm, and for the average density 
p = 4 x 10 5 g cm~ 3 Mo.6-R^ 3 - The central density is about an order of magnitude higher, 
p c = 3 x 10 6 g cm~ 3 M .6-R9 3 . Rubakov || estimated the product of cross section for catalysis 
and relative velocity v of the monopole and nucleon to be constant: av = cr = lCT 28 (cru__28) 
cm 2 . (Throughout, we take % = ks = c = 1.) For the thermal nucleon velocities expected 
inside a carbon and oxygen white dwarfs, v ~ lCT 3 c, suppression effects may reduce the 
cross section by a factor of 10~ 2 s_2 |16|], and so we include this factor. In white dwarfs made 
of helium the suppression effects would be less effective (s_2 = 10), and all the monopole 
flux bounds would be an order of magnitude stronger. 
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II. ESTIMATION OF THE MONOPOLE FLUX 



As we noted in the introduction, monopole-catalyzed nucleon decay caused by the 
monopoles captured in white dwarfs can provide an additional internal heat source for the 
star. Our evaluation of the monopole flux is based on the observed luminosity of the white 
dwarfs and estimation of the number of monopoles trapped inside the stars. As a monopole 
passes through a white dwarf, it loses energy and is captured. Electronic interactions are 
considered to be the primary source of energy loss for the monopoles, with |I7|, flgfl , 



where n e is the electron density, the Fermi momentum of the electrons kp ~ O.lMeV, 
Zmin = 2kp\/h, A is the mean free path of the electron, p is the density of the white dwarf 
(in gem -3 ), and (3 is the velocity of the monopole as it passes through the white dwarf. 

As was shown in [Q, a white dwarf accumulates almost all the monopoles with m < 10 20 
GeV incident upon it. The number of monopoles captured by a white dwarf exposed to a 
monopole flux F(cm -2 s -1 sr -1 ) for a time r = r 10 10 10 yr is given by 



where A = 4nR 2 [l + {2GM / Rv 2 M )] is the capture area, % = W-3IO 3 c is the monopole 



Once captured, the monopoles sink to the center of the white dwarf. In calculating the 
luminosity from catalyzed nucleon decay, we use the central density of the white dwarf. We 
are justified in doing this since the time for the monopole to fall (from rest) into the center is 
~ 1000s, as has been calculated |14| by treating the motion of the monopoles as a harmonic 
oscillator with a dEj dx damping term. We find the luminosity from catalyzed nucleon decay 
per monopole: 





N M « FAr(nsr) « 2.0 x 10 3 VF 



(3) 



velocity and a± = Ti R 9 M £V_l. 



Li = p c av = 8.1 x 10 7 ergs \av)^28S-2M 0S R g 3 . 



(4) 
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Then the total luminosity of a white dwarf due to a monopole-catalyzed nucleon decay 

is 



Lmon = N M Lt = 1.6 x lO^Fergs- 1 , (5) 

where a 2 = t w R^ 2 Ml_ & {av)^S-ivZl. 

From the Stefan-Boltzman law, L mon = cbb^^T^^ where obb is the Stefan-Boltzman 
constant, we can find the blackbody temperature corresponding to this luminosity, 

T eff = 1.3 x 10 8 K(Fa 2 i?9 2 )3. (6) 

White dwarfs cool as they age. The cooling time is a function of the white dwarf mass 
and composition. For white dwarf 1136-286, we use the mass and composition provided 
by the observers ||15|| . The observed energy distributions are obtained from a combination 
of both optical BVRI and infrared HJK photometric data and used to derive both the 
effective temperature and the atmospheric composition of the star. This white dwarf is 
seen to have a Helium atmosphere. Stellar masses were also obtained with trigonometric 
parallax. Bergeron, Ruiz, and Leggett |15| derive M = 1.2M for WD 1136-286. Then from 
measurements of the luminosity and T e ff, Eqn. (6) implies that the radius is R ~ 3.9xl0 8 cm. 

We also use two different cooling models. First we use the white dwarf cooling theory from 



the calculations of Segretain et al. |19 |, as communicated by G. Chabrier. The Segretain et 
al. Jl9| model accounts for gravitational energy release due to carbon-oxygen differentiation 
at crystallization. This treatment of crystallization yields significantly longer white dwarf 
cooling times, which in turn imply an older age for any particular white dwarf. These white 
dwarf models correspond to a mass sequence of initially unstratified white dwarfs composed 
of equal parts carbon and oxygen, with helium atmospheres. With these models, the age of 



white dwarf 1136-286 is 9.63 Gyr. For comparison we also use the cooling curves of Wood [20 
which do not include chemical fractionation. Chemical fractionation provides an additional 
source of energy to be radiated away; thus models that lack it cool faster. With the Wood 
cooling models, the ages of white dwarfs are somewhat younger. Hence these models give 
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younger white dwarfs that accumulate somewhat fewer monopoles and provide somewhat 
less restrictive bounds. With the Wood cooling curve, the age of white dwarf 1136-286 is 
6.47 Gyr. To illustrate the uncertainty we provide flux bounds using both possible ages, but 
note that the discrepancy is not very great. 

The cooling models discussed above do not yet have an additional heat source due to 
monopoles. If white dwarfs have indeed been accumulating monopoles, then the monopole 
contribution to the luminosity increases linearly in time, and monopole catalyzed nucleon 
decay will eventually become the dominant source of luminosity. The minimum value of the 
total luminosity must be at least as low as 10 _4 ' 94 L Q , since white dwarf WD 1136-286 with 
this luminosity has been observed to exist. Using the mass and radius discussed previously 
for this white dwarf, we then find from Equations (3-5) that 

N M < 2.2 x W 19 (av)zl 8 s~l (?) 

With the cooling curves of Segretain et al., which include the effects of chemical frac- 
tionation, the age for this particular white dwarf WD 1136-286 is given to be 9.63 Gyr as 
mentioned above. We find a flux bound 

F < 1.3 X lO^V^I^sI^gcnrVV- 1 . (8) 

With the Wood |2(J cooling curves, the age of the white dwarf is 6.47 Gyr as mentioned 
above. Then equation (6) corresponds to a flux bound 

F < 1.9 x lO^V^I^sI^cm-VV- 1 . (9) 

This bound using the Wood cooling curves is less restrictive than the one obtained using 
the Segretain et al. cooling curves. Hence, to be conservative, in Figure 1 we plot the flux 
bound of Eq. (9). Note that the monopole velocities have been determined as a function 
of monopole mass by the following equation: Vm ~3 x lCT 3 c(10 16 GeV/m) 1//2 for monopole 
mass m < 10 17 GeV and vm ~ lCT 3 c for monopoles with mass greater or equal to 10 17 GeV 
||. Thus the flux bound is flat for monopole masses greater than 10 17 GeV and drops as 
m -1 for smaller masses. This behavior can be seen in Figure 1. 
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If the monopole flux saturates the bound in equations (8) and (9), the heat release due to 
monopole-catalyzed nucleon decay would explain the dearth of white dwarfs with luminosity 
< 1O~ 5 L . That is, monopoles may be keeping white dwarfs hot. Note that the white 
dwarf luminosity due to monopole catalyzed nucleon decay scales as L mon oc t w Mq 6 . If the 
luminosity of the coolest objects we see today is in fact due primarily to the contribution from 
monopoles, then one would in principle be able to see this dependence on white dwarf mass. 
However, one would need to be able to independently measure the white dwarf luminosity, 
mass, and age in order to test this hypothesis. 



In an earlier paper Freese checked that the presence of monopoles did not drastically 
affect the properties of the white dwarf in any way. A usual white dwarf is an isothermal, 
electron degenerate object surrounded by a very thin radiative envelope. The primary 
mechanism of heat transfer through the body of the star is conduction. In the presence 
of monopoles the white dwarf remains essentially isothermal, with a radiative envelope, 
conductive main body, and convective core so that one may conclude that monopoles have 
a negligible effect on the overall structure of white dwarfs (for more details see ||14|). 

Monopole/ antimonopole annihilation: As discussed in Freese [03, monopole- 
antimonopole annihilation has no effect on the flux bound obtained in equations (8) and (9). 
There it was shown that, if the above flux bounds are satisfied, the number of monopoles 
accumulating inside the white dwarf never reaches a sufficient abundance for annihilation to 
become significant. 



III. TIGHTER BOUNDS OBTAINED IF MONOPOLES CAPTURED BY THE 
MAIN SEQUENCE PROGENITOR ARE INCLUDED: 

During its main sequence period the progenitor of the white dwarf may also have captured 
a significant number of monopoles. These additional monopoles will lead to an even tighter 
bound on the monopole flux. In order to estimate the number of monopoles captured by the 
main sequence progenitor, we must determine its mass. Unfortunately, the transformation 
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from main sequence mass to white dwarf mass is somewhat uncertain, as discussed by Wood 
|2(J. Here we use M WD = A x exp(B x M), where A x = 0.49 and B x = 0.095. For white dwarf 
1136-286 with mass Mwd = 1.2M , we find progenitor mass M = 9AM Q . 

As a monopole passes through a MS star, it loses energy. If it loses all its initial kinetic 
energy (i.e. its energy infinitely far from the star), it is captured by the star. Since the 
energy loss increases with decreasing impact parameter, the number of monopoles captured 
by a MS star exposed to a monopole flux F for a time r = tms is just the number incident 
upon the star with surface impact parameter less than some critical value, b cr it- 

Nm = (^bl lt )(nsr)[l + (^) 2 }Ft MS} (10) 

where is the monopole velocity far from the star and v esc = {2GM / R) 1 / 2 is the escape 
velocity from the star. Frieman, Freese, and Turner previously calculated numerically 
the critical impact parameter for capture. We use those results hereQ. Given the value of 
b crit , we can substitute it into the previous equation (11) to obtain the number of monopoles 
captured by the main sequence progenitor. The sum of monopoles captured by the progenitor 
plus those captured by the white dwarf cannot exceed the maximum number allowed in 
Equation (7), so that we obtain a new flux bound. Again, we have taken monopole velocities 
to be vm ~ 3 x 10 _3 c(10 16 GeV/m) 1 ' 2 for monopole mass m < 10 17 GeV and % « 10~ 3 c for 
monopoles with mass greater or equal to 10 17 GeV ||. The inclusion of monopoles captured 
by the progenitor of the white dwarf results in a bound on monopole flux that is another 



1 Recent calculations of Ahlen and Tarle [21| indicate that equation (2) for the energy loss of 
monopoles in main sequence stars must be increased by a factor of 2. Thus b cr it should be somewhat 
larger than what was calculated by Frieman, Freese, and Turner For example, for the case 
where (3 = [3 esc = 3.3 x 10~ 3 (for a 9 M Q star), rough analytic estimates indicate that the modified 
value of b CT it is roughly given by (^ i )„ ew = 0.058 + (^ 1 )oid- ^ s the calculations of [12| are 



numerical and the difference is very small, we will continue to use the prior results of [12]. In fact, 
if one were to use the newer value of b cr it, the flux bounds would be somewhat tighter. 
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order of magnitude lower. We have evaluated the new bound for both white dwarf cooling 
models. In Table 1, we have recorded, as a function of monopole mass, the value of b crit /R, 
the number of monopoles captured by the progenitor and white dwarf, and the resultant 
flux bound. The flux bound is plotted in Figure 1. 

The flux bound is most restrictive for m ~ 10 17 GeV. One can understand the reason 
for the weaker flux bounds for masses greater and less than 10 17 GeV as follows: For masses 
less than this, it is the increasing monopole velocity as a function of decreasing monopole 
mass that drives the flux bound to become weaker for smaller mass. For masses m > 10 16 
GeV, the factor [1 + (^) 2 ] in Eq. (10) is dominated by the second term (^) 2 . As the 
monopole mass decreases below 10 17 GeV, the monopole velocity increases, the term (^f^) 2 
decreases, so that in Eq. (10) the monopole flux can increase and still maintain the same 
number of monopoles in the star and hence the same luminosity. Eventually, when the mass 
drops torn < 10 16 GeV, the first term, 1, starts to dominate in the factor [1 + (^f^) 2 ], so 
that the monopole velocity becomes unimportant and the curve becomes more and more 
flat with decreasing monopole mass. This behavior can be seen in Figure 1. For monopole 
masses greater than 10 17 GeV, the flux bound also becomes weaker, this time as a function of 
increasing mass. The reason for this is as follows. These monopoles all move with the virial 
velocity of the Galaxy ~ 10 _3 c. The heavier the monopole, the harder it is to stop. Hence 
b C rit becomes smaller for heavier masses. Thus a larger monopole flux can be accomodated 
in Eq. (10) to still obtain the same number of monopoles in the star. 

IV. NEUTRON STARS WITH MAIN SEQUENCE ACCRETION: 

We also note that the dependence on monopole mass of flux bounds due to catalysis 
in neutron stars with main sequence accretion has previously been calculated incorrectly. 
In the past the bound due to catalysis in PSR 1929+10 with main sequence accretion has 
been stated as [|n]] F(crt;/10 _28 cm 2 ) < 10 -28 cm ~ 2 s _1 sr _1 . As discussed in the previous 
paragraph, the velocity dependence of monopoles of different masses determines the shape 
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of the curve of flux bounds as a function of monopole mass. As can be seen in Table 1 
and Figure 1, the correct bounds are somewhat weaker for monopole mass other than 10 17 
GeV because of the faster velocities of monopoles with smaller masses and the lower critical 
impact parameter for monopoles with larger masses. In obtaining the numbers, we have 
assumed a main sequence progenitor of 9 M & . Then the number of monopoles captured by 
the main sequence progenitor of the neutron star is the same as the number of monopoles 
captured by the white dwarf considered in this paper. 

V. CONCLUSION 

Figure 1 shows a plot of several monopole bounds: the Parker bound, the extended 
Parker bound, neutron star bounds, and the new white dwarf bound with and without main 
sequence capture. In the plots we have used the Wood cooling curves to be conservative. 
We have found that consideration of newly observed white dwarf 1136-286 with luminosity 
10 _4 ' 94 L Q and with new calculations of white dwarf cooling curves leads to a bound on the 
monopole flux that is two orders of magnitude lower than previous bounds due to catalysis 
in white dwarfs. The new bound is F(eru/lCT 28 cm 2 ) < 1.3(1.9) x lCr 20 (t>/lCr 3 c) 2 cm ~ 2 s _1 
gr -i f or Segretain [19|] (Wood [^TJ]) cooling curves respectively, where v is the monopole 
velocity. The limit is improved by including the monopoles captured by the main-sequence 
progenitor of the white dwarf: F(crt;/10~ 28 cm 2 ) < 3.5(26) x 10~ 21 cm" 2 s -1 sr _1 for 10 17 
(10 16 ) GeV monopoles with g = g^. Flux bounds for other monopole masses and parameters 
are given in Table 1. If cooler white dwarfs are discovered, a stricter bound on the monopole 
flux will result. 

We also showed that the dependence on monopole mass of flux bounds due to catalysis 
in neutron stars with main sequence accretion has previously been calculated incorrectly. 
Previously the bound due to catalysis in PSR 1929+10 with main sequence accretion has 
been stated as [|n|] .F(eru/10~ 28 cm 2 ) < 10~ 28 cm ~ 2 s _1 sr -1 . Instead, as can be seen in 
Table 1 and Figure 1, the correct bounds are somewhat weaker for monopole mass other 
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than 10 17 GeV. 

Figure Caption Bounds on the monopole flux as a function of monopole mass. The 
Parker bound [0] due to survival of the galactic magnetic field is plotted, as is the extended 
Parker bound |5| due to survival of the magnetic field early in the history of the Galaxy Mass 
density limits (Qh 2 < 1) are plotted for a uniform density of monopoles in the universe. The 
bounds due to catalysis in white dwarf WD1 136-286 as discussed in this paper are plotted; 



the plots assume the cooling curves of Wood ||20|| , and are very similar to those obtained 
using cooling curves of Segretain etal. In addition, the bounds from this white dwarf with 
main sequence accretion (WD/MS) are plotted for g = go (solid line) and for g = 2go 
(dotted line). The bounds due to calaysis in neutron star PSR 1929+10 are plotted, as are 
bounds due to this neutron star with main sequence accretion. Again the solid line is for 
g = go and the dotted line is for g = Ig^. Note that the neutron star bounds with main 
sequence accretion have dependence on the monopole mass. 
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Table 1: Bounds on the flux F of magnetic monopoles in cm _2 s _1 sr _1 . Monopoles captured 
by white dwarfs (WD) or neutron stars (NS) and their main sequence (MS) progenitors are 
included. The white dwarf cooling time is taken to be r = 9.63 Gyr. M m is the monopole 
mass in GeV, (3 is the monopole velocity, ^g 1 is the ratio of the critical impact parameter 
for a monopole in units of the radius of the main sequence star, and the monopole charge is 
g = 69e(g/<7£)) in units of the Dirac charge go- The number of monopoles captured by the 
MS progenitor and by the white dwarf are Nms and Ny/D respectively. The second to last 
column is the flux bound due to catalysis in WD 1136-286 (with MS monopoles included). 
The last column is the (corrected) flux bound due to catalysis in neutron star PSR 1929+10 
(with MS monopoles included). 
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Table 2: Same as table 1 for white dwarfs, but for cooling time r = 6.47 Gyr. 
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